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S2 ■ ABSTRACT 



We present full evolutionary calculations appropriate for the study of 
hydrogen-rich DA white dwarfs. This is done by evolving white dwarf progenitors 
t^^ . from the zero age main sequence, through the core hydrogen burning phase, the 

^ . helium burning phase and the thermally pulsing asymptotic giant branch phase 

\f\ . to the white dwarf stage. Complete evolutionary sequences are computed for a 

^ . wide range of stellar masses and for two different metallicities: Z = 0.01, which 

is representative of the solar neighborhood, and Z = 0.001, which is appropri- 
ate for the study of old stellar systems, like globular clusters. During the white 
K> , dwarf cooling stage we compute self-consistently the phase in which nuclear re- 

actions are still important, the diffusive evolution of the elements in the outer 
layers and, finally, we also take into account all the relevant energy sources in the 
deep interior of the white dwarf, like the release of latent heat and the release of 
gravitational energy due to carbon-oxygen phase separation upon crystallization. 
We also provide colors and magnitudes for these sequences, based on a new set 
of improved non-gray white dwarf model atmospheres, which include the most 
up-to-date physical inputs like the Lya quasi-molecular opacity. The calculations 
are extended down to an effective temperature of 2,500 K. Our calculations pro- 
vide a homogeneous set of evolutionary cooling tracks appropriate for mass and 
age determinations of old DA white dwarfs and for white dwarf cosmo chronology 
of the different Galactic populations. 
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Subject headings: stars: evolution — stars: white dwarfs — stars: interiors 



1. Introduction 

White dwarfs are the most common stellar evolutionary end-point. As a matter of fact, 
more than 97% of stars are expected to end their lives as white dwarfs. Their evolution 
can be described as a simple cooling process which lasts for very long periods of time. 
Moreover, since these fossil stars are abundant and long-lived objects they convey important 
information about the properties of all galactic populations (Hansen & Liebert 2003). In 
particular, white dwarfs can be considered as reliable cosmic clocks to infer the age of a 
wide variety of stellar populations, such as the Galactic disk and halo (Winget et al. 1987; 
Garcia-Berro et al. 1988a,b; Isern et al. 1998; Garcia-Berro et al. 1999; Torres et al. 2002), 
and the system of globular and open clusters (Kalirai et al. 2001; Hansen et al. 2002, 2007). 
They also provide helpful information about the past history of the star formation rate of 
our Galaxy (Diaz-Pinto et al. 1994). Additionally, they can be used to place constraints on 
exotic elementary particles (Isern et al. 1992; Corsico et al. 2001; Isern et al. 2008) or on 
alternative theories of gravitation (Garcia-Berro et al. 1995; Benvenuto et al. 2004). Recent 
reviews on the properties and evolution of white dwarfs and of their applications are those 
of Winget & Kepler (2008) and Fontaine & Brassard (2008). However, to use them for this 
wide range of applications reliable and accurate evolutionary sequences must be computed. 

Previous recent evolutionary calculations of hydrogen-rich white dwarfs are those of 
Fontaine et al. (2001), Salaris et al. (1997; 2000), Hansen (1998; 1999), Benvenuto & 
Althaus (1999) and Wood (1992; 1995). All these works have studied different aspects of the 
evolution of white dwarfs with hydrogen-rich envelopes. For instance, the most commonly 
used models, those of Wood (1995), cover a wide range of stellar masses and envelope masses 
and, until recently, were considered to be a standard reference in the field of white dwarf 
cosmo chronology. However, we emphasize that these models were computed using gray 
atmospheres, a severe drawback (especially at low luminosities) that more recent calculations 
have overcome. Among these models we mention the works of Hansen (1998; 1999). This set 
of cooling models pioneered the usage of detailed model atmospheres as surface boundary 
conditions in cooling calculations for old white dwarfs. This is an important issue since it 
affects the location of the base of the convective envelope. These calculations also showed 
that collision-induced absorption processes affect the colors of old white dwarfs. Salaris et al. 
(1997; 2000) focused on the question of the interior abundances of carbon and oxygen, which 
is of critical importance to derive reliable ages and on the effects of phase separation of carbon 
and oxygen upon crystallization, while Fontaine et al. (2001) were the first to discover the 
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importance of convective coupling between the atmosphere and the degenerate core. This last 
issue also bears importance for the determination of accurate ages. Nevertheless, all these 
works suffer from the same shortcoming. All of them evolved intial white dwarf configurations 
which were not obtained self-consistently from models evolving from the main- sequence. As 
a consequence, the chemical stratification of the inner, degenerate core was simplistic in most 
of the cases, except in the case of the cooling sequences of Salaris et al. (1997; 2000). Also, 
the envelope mass and the outer layer chemical distribution were idealized ones in all the 
cases. Additionally, the oldest sequences used physical inputs which are nowadays outdated 
and, finally, most of them disregarded the energy release of carbon-oxygen phase separation 
(Garcia-Berro et al. 1988a; 1988b). 

The aim of this paper is to compute a set of new cooling sequences for hydrogen-rich 
white dwarfs, incorporating the most up-to-date physical inputs. We emphasize that our 
evolutionary sequences are derived from a full and self-consistent treatment of the complete 
evolutionary history of progenitor stars with different masses evolved with two different 
metallicities {Z = 0.01 and Z = 0.001) appropriate for the study of the solar neighborhood 
and of metal-poor stellar systems, like globular clusters or the galactic halo. Thus, our 
calculations constitute a comprehensive set of evolutionary tracks which allow to study the 
evolution of hydrogen-rich white dwarfs in a wide variety of stellar systems. Moreover, since 
our calculations encompass the pre-white dwarf evolutionary phases, the white dwarf evolu- 
tionary calculations presented here are not affected by inconsistencies arising from artificial 
procedures to generate starting white dwarf configurations. In particular, the calculation 
of the evolutionary history of progenitor stars provides us with the amount of H left in the 
white dwarf, and with the chemical profiles expected not only for the carbon-oxygen core, 
but also for the partially degenerate regions above the core, of relevance for the white dwarf 
cooling phase. 

The paper is organized as follows. In Sect. |2] we describe the main physical inputs 
of our models. We also describe the main characteristics of the initial models, the model 
atmospheres employed in this work and some details of the evolutionary computations. In 
Sect. 3 we present the results of the pre-white dwarf evolutionary sequences, paying special 
attention to the derived initial-to-final-mass relationship, and also the white dwarf cooling 
tracks. In this section we also discuss in detail the effects of phase separation upon crystal- 
lization. Finally, in Sect. 4 we summarize our results, we discuss their significance and we 
draw our conclusions. 
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2. Computational details 

2.1. Input physics 

The evolutionary calculations presented in this work were done with an updated version 
of the LPCODE stellar evolutionary code — see Althaus et al. (2005a) and references therein. 
We used this code to compute both the evolution of white dwarfs and that of their progenitor 
stars. We emphasize that in recent years, the LPCODE stellar evolutionary code has been 
employed to study different aspects of the evolution of low-mass stars, such as the formation 
and evolution of H-deficient white dwarfs, PG 1159 and extreme horizontal branch stars 
(Althaus et al. 2005a; Miller Bertolami & Ahhaus 2006; Miller Bertolami et al. 2008; 
Althaus et al. 2009a) and, more recently, it has also been used to study the formation of hot 
DQ white dwarfs (Althaus et al. 2009b), the evolution of He-core white dwarfs with high 
metallicity progenitors (Althaus et al. 2009c), and the evolution of hydrogen-deficient white 
dwarfs (Althaus et al. 2009d). Moreover, this code has also been used to study the white 
dwarf initial-final mass relationship (Salaris et al. 2009). A recent test of the code and a 
thorough comparison of the results obtained using LPCODE with those obtained using other 
evolutionary codes has recently been made in Salaris et al. (2009). Details of LPCODE can be 
found in these works and in Althaus et al. (2009d). For these reasons, in what follows, we 
only comment on the main physical input physics, namely, those that are relevant for the 
evolutionary calculations presented in this work. 

To begin with, we recall that the radiative opacities employed in our calculations were 
those of the OPAL project (Iglesias & Rogers 1996), including carbon- and oxygen-rich 
compositions, supplemented at low temperatures with the Alexander & Ferguson (1994) 
molecular opacities. For the present calculations, we have not considered carbon-enriched 
molecular opacities (Marigo 2002), which are expected to reduce effective temperatures at the 
AGB (Weiss & Ferguson 2009). We adopted the conductive opacities of Cassisi et al. (2007), 
which cover the entire regime where electron conduction is relevant. Neutrino emission rates 
for pair, photo, and bremsstrahlung processes were those of Itoh et al. (1996), while for 
plasma processes we included the treatment presented in Haft et al. (1994). 

LPCODE considers a simultaneous treatment of non-instantaneous mixing and burning of 
elements. Specifically, abundance changes are described by the set of equations 
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where X is a vector containing the abundances of all the elements — see Althaus et al. 
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(2003) for details on the numerical procedure. In this equation, the first term describes the 
nuclear evolution, and it is fully coupled to the current composition change due to mixing 
processes, as represented by the second term. In particular, the efficiency of convective 
mixing is described by the diffusion coefficient D, which in this work is given by the mixing 
length theory. The nuclear network accounts explicitly for the following 16 elements: ^H, 
2H, 3He, ^He, ^Li, ^Be, ^^C, i^c, i^n, i5n, isq, itq, i^O, ^^F, ^ONe and 22Ne, together with 
34 thermonuclear reaction rates of the pp chains, CNO bi-cycle, helium burning, and carbon 
ignition that are identical to those described in Althaus et al. (2005a), with the exception of 
the reaction ^^C + p — )■ ^^N + 7 — )• ^^C + e^ + v^ and ^^C(p, 7)^'^N, for which we adopted the 
rate of Angulo et al. (1999). In particular, it is worth noting that the ^^C(a;,7)^^0 reaction 
rate — which is of special relevance for the carbon-oxygen stratification of the emerging white 
dwarf core — adopted here was also that of Angulo et al. (1999). In passing, we mention 
that a detailed inner composition is required for a proper computation of the energy released 
by the redistribution of chemical elements during crystallization of the white dwarf core (see 
below). Specifically, the energy released during this process, and the resulting impact on 
the cooling ages of faint white dwarfs, increases for larger carbon abundances (Isern et al. 
1997; 2000). The standard mixing length theory for convection — with the free parameter 
a = 1.61 — was adopted. With this value, the present luminosity and effective temperature 
of the Sun, log Teg = 3.7641 and Lq = 3.842 x 10^^ erg s-\ at an age of 4570 Myr, are 
reproduced by LPCODE when Z = 0.0164 and X = 0.714 are adopted — in agreement with 
the Z/X value of Grevesse & Sauval (1998). 

During the evolutionary stages prior to the thermally pulsing AGB (TP-AGB) phase, we 
allowed the occurrence of extra-mixing episodes beyond each convective boundary following 
the prescription of Herwig et al. (1997). As well known, the occurrence of extra-mixing 
episodes is of relevance for the final chemical stratification of the white dwarfs, particularly 
during the late stage of core helium burning phase — see Prada Moroni & Straniero (2002) 
and Straniero et al. (2003). We treated extra- mixing as a time- dependent diffusion process 

— by assuming that mixing velocities decay exponentially beyond each convective boundary 

— with a diffusion coefficient given by Z^em = -Do Gxp{—2z/fHp), where Hp is the pressure 
scale height at the convective boundary, Dq is the diffusion coefficient of unstable regions 
close to the convective boundary, and z is the geometric distance from the edge of the 
convective boundary (Herwig et al. 1997). For all our sequences we adopted / = 0.016, a 
value inferred from the width of the upper main sequence. In the case of the sequence of 
1.25 Mq with Z = 0.001, the size of the convective core on the main sequence is small. In 
this case we use / = 0.008 during the core hydrogen burning phase. Abundance changes 
resulting from extra- mixing were fully coupled to nuclear evolution, following Eq. ([T]). 

In the present study, extra-mixing episodes were disregarded during the TP-AGB phase. 
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In particular, a strong reduction (a value of / much smaller than 0.016) of extra-mixing 
episodes at the base of the pulse-driven convection zone seems to be supported by simu- 
lations of the s— process abundance patterns (Lugaro et al. 2003) and, more recently, by 
observational inferences of the initial-final mass relation (Salaris et al. 2009). As a result, it 
is expected that the mass of the hydrogen-free core of our sequences gradually grows as evo- 
lution proceeds through the TP-AGB. This is because a strong reduction or suppression of 
extra-mixing at the base of the pulse-driven convection zone strongly inhibits the occurrence 
of third dredge-up, thus favoring the growth of the hydrogen-free core. The implications 
of this treatment for the theoretical initial-final mass relation will be discussed later in this 
paper — see also Salaris et al. (2009) and Weiss & Ferguson (2009). The breathing pulse 
instability occurring towards the end of core helium burning was suppressed — see Straniero 
et al. (2003) for a discussion of this point. 

Mass loss was considered during core helium burning and red giant branch phases fol- 
lowing Schroder & Cuntz (2005). During the AGB and TP-AGB phases, we considered the 
prescription of Vassiliadis & Wood (1993). In the case of a strong reduction of third dredge- 
up, as occurred in our sequences, mass loss plays a major role in determining the final mass 
of the hydrogen-free core at the end of the TP-AGB evolution, and thus the initial-final mass 
relation. In particular, mass loss erodes the hydrogen-rich envelope of the star and limits 
the additional growth of the core mass during the TP-AGB. This is quite different from the 
situation in which appreciable third dredge-up takes place, in which case the final core mass 
at the TP-AGB phase is not very different from the mass at the first thermal pulse (Weiss 
& Ferguson 2009), and the role of mass loss becomes less relevant. 

We want to mention that LPCODE is able to follow the complete evolution of the star 
model from the main sequence through the core helium fiash and, finally, for many thermal 
pulses on the AGB almost without hand intervention, except for the latest stages of the 
TP-AGB phase of the more massive sequences, where numerical instabilities arise. These 
instabilities, also found by other authors, are the result of the large radiation pressure in the 
convective envelope. To circumvent them, we found computationally convenient to artificially 
modify the opacity profile in those regions close to the base of the convective envelope. This 
procedure bears no relevance for the evolution during the white dwarf regime. 

For the white dwarf regime, the physical inputs considered in LPCODE were completely 
revised and updated from those considered in our previous studies, particularly regarding the 
treatment of crystallization. With regard to the microphysics, for the high-density regime, 
we used the equation of state of Segretain et al. (1994), which accounts for all the important 
contributions for both the liquid and solid phases — see Althaus et al. (2007) and references 
therein. For the low-density regime, we used an updated version of the equation of state 
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of Magni & Mazzitelli (1979). Radiative and conductive opacities are those described at 
the beginning of this section. In particular, conductive opacities are considered for densities 
larger than that at which the main chemical constituents are completely pressure-ionized. 
During the white dwarf regime, the metal mass fraction Z in the envelope is not assumed to 
be fixed, rather, it is specified consistently acording to the prediction of element diffusion. To 
account for this, we considered radiative opacities tables from OPAL for arbitrary metallici- 
ties. For effective temperatures less than 10,000K we include the effects of molecular opacity 
by assuming pure hydrogen composition from the computations of Marigo & Aringer (2009). 
This assumption is justified because element diffusion leads to pure hydrogen envelopes in 
cool white dwarfs. For the white dwarf regime, convection was treated in the formalism of 
the mixing length theory, as given by the ML2 parameterization (Tassoul et al. 1990). 

As well known, there are several physical processes that change the chemical abundance 
distribution of white dwarfs along their evolution. In particular, element diffusion strongly 
modifies the chemical composition profile throughout their outer layers. In this work, we 
computed the white dwarf evolution in a consistent way with the changes of chemical abun- 
dance distribution caused by element diffusion along the entire cooling phase. As a result, 
our sequences developed pure hydrogen envelopes, the thickness of which gradually increases 
as evolution proceeds. We considered gravitational settling and thermal and chemical diffu- 
sion — but not radiative levitation, which is only relevant at high effective temperatures for 
determining the surface composition — of ^H, ^He, *^He, ^^C, ^^C, ^^N and ^^O, see Althaus 
et al. (2003) for details. In particular, our treatment of time-dependent diffusion is based 
on the multicomponent gas treatment presented in Burgers (1969). In LPCODE, diffusion 
becomes operative once the wind limit is reached at high effective temperatures (Unglaub & 
Bues 2000). We assume this to occur when the surface gravity of our models g > gdis, where 
^diff = 7 if Tcff > 90, OOOK and c/dis = 6.4 + Teg/ISO, OOOK if T^s < 90, OOOK. For smaller 
gravities, wind mass-loss is high enough that prevents appreciable element diffusion from 
occurring. This prescription represents reasonably well the detailed simulations of Unglaub 
& Bues (2000) for the occurrence of wind limits in hydrogen-rich white dwarfs. Other phys- 
ical process responsible for changes in the chemical profile of white dwarfs that we took into 
account is related to carbon-oxygen phase separation during crystallization. In this case, 
the resulting release of gravitational energy considerably impacts the white dwarf cooling 
times. Abundance changes resulting from residual nuclear burning — mostly during the 
hot stages of white dwarf evolution — and convective mixing, were also taken into account 
in our simulations. In particular, the release of energy by proton burning was considered 
down to log(L/LQ) ^ —4. The role of residual hydrogen burning in evolving white dwarfs 
is by no means negligible, particularly in the case of those white dwarfs resulting from low- 
metallicity progenitors. Finally, we considered the chemical rehomogenization of the inner 



carbon-oxygen profile induced by Rayleigh- Taylor instabilities following Salaris et al. (1997). 
These instabilities arise because the positive molecular weight gradients that remain above 
the fiat chemical profile left by convection during helium core burning. 

An important aspect of the present work was the inclusion of energy sources resulting 
from the crystallization of the white dwarf core. This comprises the release of latent heat 
and the release of gravitational energy associated with changes in the carbon-oxygen profile 
induced by crystallization. In this study, the inclusion of these two energy sources was done 
self-consistently and locally coupled to the full set of equations of stellar evolution. That 
is, we computed the structure and evolution of white dwarfs with the changing composition 
profile and with the luminosity equation appropriately modified to account for both the 
local contribution of energy released from the core chemical redistribution and latent heat. 
This constitutes an improvement over previous attempts (Salaris et al. 2000) to include 
the release of energy from phase separation in stellar evolutionary codes. Details about the 
numerical procedure to compute the energy sources from crystallization will be presented in 
a forthcoming work. Briefly, at each evolutionary timestep we computed the crystallization 
temperature and the change of chemical composition resulting from phase separation using 
the spindle-type phase diagram of Segretain & Chabrier (1993). This phase diagram provides 
the crystallization temperature (which depends on the chemical composition of the liquid 
phase) as a function of the crystallization temperature of a one component plasma. In 
our calculations, the one component plasma crystallization temperature is computed by 
imposing T = 180, where T = {Z^/'^)e'^ /a^k-oT is the ion coupling constant, and ag is the 
interelectronic distance. After computing the chemical compositions of both the solid and 
liquid phases we evaluated the net energy released in the process as in Isern et al. (1997), and 
added to it the latent heat contribution, of the order of 0. 77 k-sT per ion, which is usually 
smaller. Both energy contributions were distributed over a small mass range around the 
crystallization front. We mention that the magnitude of both energy sources was calculated 
at each iteration during the convergence of the model. 



2.2. Model atmospheres 

A proper treatment of the evolutionary behavior of cool white dwarfs requires the use 
of outer boundary conditions as provided by detailed non-gray model atmospheres. To this 
end, we considered model atmospheres that incorporate non-ideal effects in the gas equation 
of state and chemical equilibrium, collision-induced absorption from molecules, and the Lya 
quasi-molecular opacity. Specifically, for Tcs < 10, OOOi^, we derived starting values of the 
pressure, temperature, radial thickness and outer mass fraction at an optical depth r = 25 
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from a grid of non-gray model atmospheres which covers a surface gravity range between 
log^f = 6.5 and 9. For larger values of r, the use of Rosseland mean opacities is justified, 
and the diffusion approximation for the radiative transfer can be assumed. The use of 
non-gray model atmospheres to derive outer boundary conditions gives rise to shallower 
outer convection zones, as compared with the standard gray treatment of the atmosphere 
(Bergeron et al. 1997). At advanced stages of white dwarf evolution, the central temperature 
becomes strongly tied to the temperature stratification of the outer layers. Thus, using non- 
gray model atmospheres is highly desired for an accurate assessment of cooling times of cool 
white dwarfs (Prada Moroni & Straniero 2007). 

In the present work, model atmospheres were specifically computed on the basis of im- 
proved LTE model atmospheres. Colors and magnitudes were evaluated for effective temper- 
atures lower than 60,000K, because NLTE effects become important above this temperature. 
Calculations were done for a pure hydrogen composition and for the HST ACS filters (Vega- 
mag system) and UBVRI photometry. The numerical code used is a new and updated 
version of the one described in Rohrmann et al. (2002). Models were computed assuming 
hydrostatic and radiative-convective equilibrium. Convective transport present in the cooler 
atmospheres was treated within the usual mixing-length (ML2) approximation, adopting 
the same value of a used in the evolutionary calculations. The microphysics included in 
the model atmospheres comprises non-ideal effects in the gas equation of state and chemical 
equilibrium based on the occupation probability formalism as described in Rohrmann et al. 
(2002). The code includes H, H2, H+, H", H+, H+, He, He", He+, He2+, He^, HeH+, and 
e~. The level occupation probabilities are self-consistently incorporated in the calculation of 
the line and continuum opacities. Collision-induced absorptions due to H2-H2, H2-He, and 
H-He pairs are also taken into account (Rohrmann et al. 2002). 

For the purpose of the present work, the model atmospheres explicitly included the 
Lya quasi-molecular opacity according to the approximation used by Kowalski & Saumon 
(2006). Quasi-molecular absorption results from perturbations of hydrogen atoms by in- 
teractions with other particles, mainly H and H2. Here, we considered extreme pressure- 
broadening of the line transition H(n = 1)— 7'H(n = 2) due to H-H and H-H2 collisions, with 
the red wing extending far into the optical region. A detailed description of these collisional 
line-broadening evaluations will be presented in a forthcoming paper. On the basis of the 
approximations outlined in Kowalski & Saumon (2006), we evaluated the red wing absorp- 
tion within the quasi-static approach using theoretical molecular potentials to describe the 
interaction between the radiator and the perturber. We also considered the variation in the 
electric-dipole transition moment with the interparticle distance. The H3 energy-potential 
surfaces contributing to collisions H-H2 were taken from Kulander & Guest (1979) and Roach 
& Kuntz (1986), and the dipole transition moments were calculated from Petsalakis et al. 
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(1988). Broadening of Lya line by H-H collisions plays a minor role compared to H-H2 
encounters. The potential interactions for H-H configurations were taken from Kolos & Wol- 
niewicz (1965) and the transition probability was assumed constant in this case. The main 
effect of the Lya quasi-molecular opacity is a reduction of the predicted flux at wavelength 
smaller than 4000 A for white dwarfs cooler than Teg ~ 6, 000 K. 



2.3. Initial models 

The initial models for our white dwarf sequences correspond to stellar configurations 
derived from the full evolutionary calculations of their progenitor stars. The initial He 
content of our starting models at the main sequence were provided by the relation Y = 
0.23 + 2.41Z, as given by present determinations of the chemical evolution of the Galaxy 
(Flynn 2004; Casagrande et al. 2007). Two metallicities for the progenitor stars were 
considered: Z = 0.01 and 0.001. Hence, the initial compositions of our sequences are, 
respectively, (Y, Z) = (0.254,0.01) and {¥, Z) = (0.232,0.001). All the sequences were 
evolved from the ZAMS through the thermally-pulsing and mass-loss phases on the AGB 
and, finally, to the white dwarf cooling phase. Specifically, we computed 10 full evolutionary 
sequences for Z = 0.01 and six for Z = 0.001. In Table [H we list the initial masses of 
the progenitor stars at the ZAMS, together with other evolutionary quantities which will be 
discussed below. 



3. Evolutionary results 

3.1. Prom the ZAMS to the white dwarf stage 

The evolution in the Hertzsprung-Russell diagram of our sequences from the ZAMS to 
advances stages of white dwarf evolution is shown in Fig. [H for the case in which Z = 0.01 
is adopted. Note that the less massive sequence experiences a hydrogen subflash before 
entering its final cooling track. The initial masses at the ZAMS and the final white dwarf 
masses of these sequences can be found in Table [H for both metallicities. In this table we 
also list the main-sequence lifetimes, which, as well known, for solar-metallicity sequences 
are larger than those of their metal-poor counterparts. We mention that our models have 
main-sequence lifetimes longer than those recently published by Weiss & Ferguson (2009). 
Differences are less than 8 % however, except for the main-sequence lifetime of the solar 
sequence, which is 25 % larger. One of the reason for such a discrepancy is due to our 
simplified treatment of the equation of state during the evolutionary stages prior to white 
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dwarf formation. Also listed in this table are the total masses of the residual hydrogen 
and helium content left in the white dwarfs at the evolutionary stage corresponding to the 
point of maximum effective temperature in the Hertzsprung-Russell diagram. Note that the 
residual hydrogen content decreases with the white dwarf mass, a well-known result. For 
the case of solar metallicity progenitors, the hydrogen mass differs by a factor of 20 for 
the stellar mass range considered. This general trend is also observed for the mass of the 
helium content, where the mass of the residual helium ranges from 0.025 to 0.0022 Mq. The 
hydrogen and helium masses listed in Table [T] should be considered as upper limits for the 
maximum hydrogen and helium content left in a white dwarf resulting from the evolution 
of single star progenitors. However, the occurrence of a late thermal pulse after departure 
from the TP-AGB may reduce the hydrogen mass considerably, see Althaus et al. (2005b). 

Note also that for a given white dwarf mass there is a marked dependence of the final 
hydrogen mass on the initial metallicity of the progenitor star: higher hydrogen masses 
are expected in metal-poor progenitors, see Iben & MacDonald (1985, 1986). For instance, 
for the 0.593 M0 white dwarf model we find that the hydrogen mass is logMn — —3.950 
for the solar metallicity progenitor, while this mass turns out to be logMn — —3.777, for 
the metal-poor progenitor, i.e. a 50% higher. This is an important issue since one of the 
factors affecting the white dwarf cooling rate is, precisely, the thickness of the hydrogen-rich 
envelope. By contrast, in the case of the residual helium content, no appreciable dependence 
on the metallicity exists. Also shown in Tabled] are the number of thermal pulses during the 
AGE and the total mass lost during the entire AGB phase, in solar units. We find that the 
number of thermal pulses during the AGB phase is, generally speaking, slightly larger for 
the set of metal-poor evolutionary sequences. For instance, for the 1.5 M0 stellar sequence 
a total of 7 thermal pulses occur for the solar metallicity white dwarf progenitor, while for 
the metal-poor progenitor this number is 10, thus leading to a more extended mass-loss 
phase. However, the total mass lost during the entire AGB phase is smaller for the case of 
a metal-poor progenitor — 0.757 Mq for the same model star — than for a solar-metallicity 
progenitor — 0.795 Mq. 

Perhaps one of the most interesting results of our full evolutionary calculations is the 
initial-to-final mass relationship. In Fig. H] we show our results for the case of solar com- 
position {Z = 0.01). Specifically, in this figure we show using a solid line the mass of the 
white dwarf resulting from our theoretical calculations as a function of the initial mass at the 
ZAMS. For the sake of comparison we also show using a dot-dashed-dashed line the mass 
of the hydrogen-free core at the beginning of the first thermal pulse as a function of the 
initial mass of the progenitor star, and we compare this relationship with that of Salaris et 
al. (1997) — short dashed line — and Dominguez et al. (1999) for Z = 0.006 — long dashed 
line. These two initial-to-final mass relationships were obtained assuming that the mass of 
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the resulting white dwarf corresponds to that of the hydrogen-free degenerate core at the 
first thermal pulse, and consequently that the core does not grow appreciably afterwards. 
In addition, we also show (using a dot-dot-dashed line) the recently obtained initial-to-final 
mass relationship of Weiss & Ferguson (2009) for Z = 0.008. In this calculation envelope 
overshooting was used during the TP-AGB, which, as mentioned in Sect. 12.11 considerably 
reduces the further growth of the hydrogen-free core. We also show the observational initial- 
to-final mass relationship of Catalan et al. (2008a) — dot-dashed line. This relationship is 
based on cluster observations of different metallicities, which are close to that adopted in 
our calculations {Z = 0.01). Also shown are the observational results of Kalirai et al. (2008) 
— with their corresponding error bars — which correspond to metallicities very close to the 
solar metallicity — filled and open circles — and the results for individual white dwarfs in 
common proper motion pairs of Catalan et al. (2008b) — solid triangles. The main sequence 
stars of these common proper motion pairs also have solar metallicity. In addition, the re- 
sults for M35 as quoted in Catalan et al. (2008a), for which the estimated metallicity is also 
close to solar (Barrado y Navascues et al. 2001) are also shown. Finally, the semi-empirical 
initial-to-final mass relationship of Safaris et al. (2009) based on open cluster observations 
is included in Fig. HJ 

It is to be noted the excellent agreement between our theoretical calculations and the 
empirical initial-to-final mass relationships, particularly that of Salaris et al. (2009). Note 
as well that the mass of the hydrogen-free degenerate core at the first thermal pulse for all 
the theoretical sequences agrees with each other and presents a minimum around ~ 2.0 Mq, 
but does not agree with the empirical initial-to-final mass relationship. This emphasizes the 
importance of carefully following the evolution of the star models from the main sequence 
all the way through the TP-AGB phase and, finally, to the beginning of the white dwarf 
cooling track, when the mass-loss rate becomes negligible. In particular, the growth of the 
core mass during the TP-AGB phase is emphasized as a gray area in Fig. [2J The implica- 
tion of a proper computation of the intial-to-final mass relationship for the carbon/oxygen 
composition expected in a white dwarf will be discussed in a forthcoming paper. Note as 
well, that our pre-white dwarf evolutionary calculations provide us with accurate and reliable 
starting configurations at the beginning of the white dwarf cooling phase, as they yield not 
only self-consistent inner chemical profiles, but also masses of the hydrogen-rich envelopes, 
helium buffers and core masses which agree with the observational results. 
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3.2. A global view of the white dwarf cooling phase 

In Fig. |3]we show the different luminosity contributions during the white dwarf cooling 
phase, for an archetypical 0.609 Mq carbon-oxygen white dwarf resulting from a progeni- 
tor star of 2.0 Mq with solar composition. The very first phases of the cooling phase are 
dominated by residual hydrogen burning in the outer layers. This can be easily seen in 
Fig. 131 where the different nuclear luminosities, namely the proton-proton hydrogen-burning 
luminosity, the CNO bicycle luminosity and the helium-burning luminosity are shown as a 
function of the cooling age. As can be seen in this figure, at high luminosities the largest con- 
tribution comes from the CNO bicycle, being the proton-proton and the helium-burning lu- 
minosities orders of magnitude smaller. This a short-lived phase (a few thousand years) and, 
thus, given the long-lived cooling times of white dwarfs, it is totally negligible in terms of age. 
Nevertheless, this phase is important as it configures the final thickness of the hydrogen-rich 
envelope of the white dwarf. After this short-lived phase, the nuclear luminosities abruptly 
decline (at logt ~ 3.6) and the release of gravothermal energy becomes the dominant energy 
source and drives the evolution. In passing, we note that residual nuclear reactions are not 
totally extinguished until very late phases of the evolution, in agreement with the pioneer- 
ing results of Iben & MacDonald (1985, 1986). In fact, there are still small contributions 
of both the CNO cycle and proton-proton chains until logt ~ 8.3 and logt ~ 9.0, respec- 
tively. Although almost negligible for the calculation of the cooling age in the case of white 
dwarfs resulting from solar metallicity progenitors, this residual nuclear burning becomes 
relevant at very late stages for white dwarfs resulting from low-metallicity progenitors as it 
lasts for one billion years, see next section. The phase in which the evolution is driven by 
gravitational contraction lasts for about one million years. During this phase the release of 
gravothermal energy occurs preferentially in the outer partially- or non-degenerate layers of 
the white dwarf. More or less at the same epoch — that is, at log(t) ~ 5.6 — neutrino 
losses become also important. In particular, at this epoch neutrinos are the dominant en- 
ergy sink in the degenerate core of the white dwarf, and their associated luminosity becomes 
larger than the optical luminosity. In fact, during a relatively long period of time (from 
logt ~ 5.6 to 7.1) the neutrino luminosity is comparable to the luminosity associated to the 
gravothermal energy release. It is also interesting to note that at approximately the same 
time, element diffusion is operating in the outer partially-degenerate envelope, shaping the 
chemical stratification of the very outer layers of the white dwarf. We will discuss below the 
resulting chemical stratification. This phase lasts for about 2.2 x 10® years. At logt ~ 7.8 the 
temperatures in the degenerate core decrease below the threshold where neutrino emission 
ceases and, consequently, the neutrino luminosity abruptly drops. During this phase of the 
evolution most of the energy released by the white dwarf has gravothermal origin, and the 
white dwarf cools according to the classical Mestel's law (Mestel 1952). Finally, at logt ~ 9.2 
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crystallization sets in at the center of the white dwarf and the coohng process slows down 
due to the release of latent heat and of gravitational energy due to carbon-oxygen phase 
separation. These physical processes are noticeable as a change in the slope of the cooling 
curve. Note as well that during the crystallization phase the surface energy is larger than 
the gravothermal luminosity, a consequence of these two additional energy sources. This 
phase lasts for ~ 9.4 x 10^ years. After this phase, the temperature of the crystallized core 
of the white dwarf drops below the Debye temperature and the heat capacity of the white 
dwarf drops. Consequently, the white dwarfs enters the so-called Debye cooling phase, and 
the slope of the cooling curve increases again. This occurs at logt ~ 10. 



3.3. The thickness of the hydrogen envelope 

Fig. m shows the temporal evolution of the masses of the hydrogen content for two 
representative white dwarf cooling tracks of the two metallicities explored here. Also shown 
are the ratio of the hydrogen-burning luminosities to the total luminosity. In particular, 
the thick lines represent the evolution of a 0.609 M© white dwarf resulting from a solar 
progenitor, while the thin lines show that of a 0.593 Mq white dwarf resulting from a metal- 
poor progenitor. The solid lines correspond to the evolution of the mass of the hydrogen 
content, while the dashed lines show the evolution of the nuclear luminosities. As can be 
seen, residual hydrogen burning is dominant during the first evolutionary phases of the white 
dwarf stage. As a consequence, the mass of these envelopes decreases for a period of time of 
~ 3 X 10^ years, during which hydrogen burning supplies most of the surface luminosity of the 
white dwarf. However, as soon as the mass of the hydrogen content decreases below a certain 
threshold (~ 8 x 10^^ M0 for the solar-metallicity star and ~ 1.7 x 10^^ M© for the metal- 
poor star) the pressure at the bottom of the envelope is not large enough to support further 
nuclear reactions, and hence the main energy source of the white dwarf is no longer nuclear 
reactions but gravothermal energy release, and the hydrogen content reaches a stationary 
value. This is true in the case of the white dwarf with a solar metallicity progenitor, but for 
the white dwarf remnant that results from the lower metallicity progenitor star, it is clear that 
residual hydrogen burning is by no means negligible. In this case, note that hydrogen burning 
represents an importante fraction of the surface luminosity after k. 10^ yr of evolution, and 
even at more advanced stages (~ 10^ yr), this contribution reaches up to 30 %. At this 
time, the nuclear energy production is almost entirely from the proton-proton chain, and the 
hydrogen content has been reduced down to ~ 1.1 x 10"'^ Mq. The contribution of hydrogen 
burning to surface luminosity increases for white dwarfs with lower stellar masses. We would 
like to emphasize at this point the importance of computing self-consistently the very first 
stages of the white dwarf evolution, as they provide an homogeneous set of white dwarf 
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envelope masses, which as the evolution proceeds influence the cooling of white dwarfs. 



3.4. The chemical abundances of the envelope 

Fig. [5] shows the chemical profiles of the 0.609 Mq white dwarf resulting from a solar 
metallicity progenitor for selected evolutionary stages along the white dwarf cooling track. 
Each of the panels is labelled with the luminosity and effective temperature of the evolu- 
tionary stage. In these panels we show the abundance profiles of hydrogen, helium, carbon, 
nitrogen and oxygen in terms of the outer mass fraction. As can be seen in the upper-left 
panel, which depicts the chemical profiles at the beginning of the cooling track, the resulting 
white dwarf has a hydrogen-rich envelope, with substantial amounts of heavier elements, like 
helium, carbon, nitrogen and oxygen. The chemical composition of this layer is similar to 
that of typical AGB stars that have not experienced third dredge-up episodes, being oxygen 
more abundant than carbon, and nitrogen almost as abundant as carbon. Specifically, these 
abundances are essentially fixed by the first dredge-up episode during the red giant phase. 
The deeper layers in the helium buffer zone show CNO abundances that reflect the occur- 
rence of hydrogen burning in prior stages, with nitrogen far more abundant than carbon and 
oxygen. As the white dwarf evolves across the knee in the Hertzsprung-Russell diagram, 
gravitational settling and diffusion become the relevant physical processes and the heavier 
chemical elements begin to sink appreciably. This is illustrated in the upper-right panel 
of Fig. |5l As can be seen in this panel, at this stage, the white dwarf has already devel- 
oped a thin pure hydrogen envelope that thickens as evolution proceeds. Note that at this 
evolutionary stage some diffusion and gravitational settling has already occurred in deeper 
layers, and the chemical interfaces exhibit less sharp discontinuities. During these stages, 
the chemical composition has also changed as a result of nuclear burning via the CN cycle at 
the base of the hydrogen envelope. With further cooling — see the bottom-left panel of Fig. 
|5] — the action of element diffusion becomes more apparent. In fact, the helium-rich buffer 
increases its size and both carbon and oxygen sink towards deeper and deeper regions of the 
white dwarf. Also, the thickness of the hydrogen rich layer increases appreciably, and at the 
same time, the tail of the hydrogen distribution continues to chemically diffusing inward. At 
this stage, which correponds to the domain of the pulsating DA white dwarfs, a rather thick 
hydrogen envelope has been formed, and below it, a helium-rich and several very thin layers, 
which are rich in even heavier elements — a consequence of the high gravity of the white 
dwarf. Finally, the bottom-right panel depicts the situacion after the onset of crystallization. 
Note the change of carbon and oxygen composition of the core as a result of crystallization. 
This sequence of figures emphasizes the importance of a proper treatment of time- dependent 
diffusion processes during white dwarf evolution, and the extent to which the initial chemical 
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stratification at the start of the coohng phase is ahered by these processes. 



3.5. Convective coupling and crystallization 

As discussed in Sect. 3.2, the cooling curve is influenced strongly by crystallization. 
However, at this evolutionary stage the slope of the cooling curve is not only dictated by the 
release of latent heat and other energy sources associated to crystallization but also by the 
so-called convective coupling. When the envelope becomes fully convective the inner edge of 
the convective region reaches the boundary of the degenerate regions (Fontaine et al. 2001). 
This effect is illustrated in Fig. El where we show as a function of the surface luminosity 
(in solar units) the evolution of the cooling times and central temperatures (left scales), 
and the mass of the crystallized white dwarf core (right scale) for two white dwarfs, a low- 
mass white dwarf of 0.525 M© and a rather massive white dwarf with M = 0.878 Mq, both 
resulting from solar metallicity progenitors. As can be seen there, as both white dwarfs cool, 
there is a gradual decrease of the central temperatures, while their corresponding cooling 
times also increase smoothly. At approximately log(L/L0) c^ —3.8 crystallization sets in for 
the less massive white dwarf, whereas for the massive star this occurs at log(L/L0) c^ —3. 
Convective coupling between the degenerate core and the partially degenerate convective 
envelope also occurs at low luminosities. Since the inner edge of the convective envelope 
reaches the boundary of the core, an increase of the rate of energy transfer across the outer 
opaque envelope occurs, which is much more efficient than radiative transfer alone. As 
a consequence, the relation between the central temperature and the surface luminosity 
experiences a sudden change of slope, which can be clearly seen in Fig. |6l where we show 
the region in which convective coupling occurs as a shaded area. Note that in the case of the 
more massive white dwarf, convective coupling takes place at luminosities markedly lower 
than that at which crystallization starts in the core. In fact, more than 90 % of the white 
dwarf mass has crystallized by the time convective coupling occurs in the 0.878 M© white 
dwarf. In contrast, for the less massive white dwarf, both crystallization and convective 
coupling occur at approximately the same stellar luminosity, and thus the resulting impact 
of these effects on the rate of cooling is more noticeable in this case. 



3.6. Cooling times and chemical composition of the core 

One of the most noticeable features of the white dwarf cooling tracks presented here 
is the inclusion in a self-consistent manner of the release of gravitational energy due to 
phase separation of carbon and oxygen upon crystallization. Previous studies of this kind 
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(Salaris et al. 2000) included the effects of phase separation, but a semianalytical approach 
was used. To highhght the importance of phase separation upon crystalhzation we have 
computed two different sets of coohng sequences. In the first of these coohng sequences phase 
separation of carbon and oxygen was fully taken into account, whereas in the second it was 
disregarded. In Table [2] we list for various luminosities the cooling ages of all our white dwarf 
sequences resulting from solar metallicity progenitors, when carbon-oxygen phase separation 
is neglected (top section), and the corresponding delays introduced by carbon-oxygen phase 
separation. We also show the same quantities for the case of metal-poor progenitors in 
Table |3l Clearly, phase separation of carbon and oxygen introduces significant delays at low 
luminosities, between 1.0 and 1.8 Gyr. It is worth mentioning that at log(L/LQ) = —4.6, 
a luminosity slightly smaller than that of the observed drop-off in the disk white dwarf 
luminosity function, log(L/LQ) ~ —4.5, phase separation of carbon and oxygen represents 
a correction of ~ 15% to the total age, that although not very large it is not negligible 
whatsoever if precise cosmochronology is to be performed. At log(L/L0) = —4.0, the delays 
constitute 20 — 25% of the age for the more massive white dwarfs. Note that the magnitude 
of the delays increases with the mass of the white dwarf. For the case of white dwarfs 
resulting from metal poor progenitors (see Table [3]), and for the same fiducial luminosity 
the delays introduced by carbon-oxygen phase separation are slightly larger for the same 
white dwarf mass. Our computed delays are larger than those obtained by Salaris et al. 
(2000). For instance, for our 0.609 M© white dwarf cooling sequence at log(L/L0) = —4.6 
we obtain 5t ~ 1.38 Gyr, while Salaris et al. (2000) at the same luminosity obtain for 
their 0.61 Mq white dwarf evolutionary sequence 5t ~ 1.00 Gyr. This difference stems in 
part from the larger carbon abundances of our white dwarf model, which leads to a larger 
energy release of the carbon-oxygen phase separation process, and consequently to larger time 
delays. Indeed, the chemical profiles used by Salaris et al. (2000) were those of Salaris et al. 
(1997). The central carbon abundance for the 0.61 M© white dwarf is Xq ~ 0.25, while for 
our 0.609 Mq model we obtain Xq ^ 0.29. Hence, the delays introduced by carbon-oxygen 
phase separation are correspondingly larger in our model. A realistic core composition is 
crucial for a proper assessment of the energy release from phase separation and its impact 
on the cooling times. 

Fig. [7] shows the evolutionary cooling sequences of several selected white dwarfs. Specif- 
ically we show the luminosity as a function of the cooling age for white dwarfs with so- 
lar metallicity progenitors and masses 0.525 Mq (upper-left panel), 0.570 M0 (upper-right 
panel), 0.609 Mq (bottom-left panel), and 0.877 M© (bottom-right panel), respectively. The 
figure emphasizes the evolutionary stages where the processes of convective coupling, cyrstal- 
lization and Debye cooling take place. As mentioned, in low-mass white dwarfs, cyrstalliza- 
tion and convective coupling occur approximately at the same luminosity, log(L/L0) ^ —4, 
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thus resulting in a pronounced impact on the cooling rate. As can be observed, the cooling 
tracks presented here have been computed down to very low luminosities, typically of the 
order of 10~^ Lq, or even smaller. At these very low luminosities the central regions of the 
white dwarf have low enough temperatures to enter the so-called Debye cooling phase. In 
this phase the specific heat drops abruptly as T^, and as a consequence the cooling rate is 
enhanced. Thus, the cooling curve rapidly drops. The transition to this phase of cooling 
depends on the Debye temperature, Od, which in turn depends on the density (6*0 oc p^^^). 
Consequently, more massive white dwarfs enter this phase at larger temperatures (luminosi- 
ties). In our most massive sequence, rapid Debye cooling is expected to occur at the lowest 
luminosities we computed, as it is clear from Fig. [71 while for lower stellar masses, this phase 
is delayed to much lower luminosities. 

Phase separation in the deep interiors of white dwarfs also has obvious imprints in the 
chemical profiles of carbon and oxygen. To illustrate this, in Fig. |8] we display for four 
selected white dwarf evolutionary sequences the oxygen mass abundance as a function of 
the interior mass at three selected evolutionary stages. Specifically, we show the abundance 
profiles shortly after the progenitor star departs from the AGB (dashed lines), the same 
profiles after Rayleigh- Taylor rehomogeneization (Salaris et al. 1997) has occurred (solid 
lines), and finally, when the entire white dwarf core has crystallized (dot-dot-dashed lines). 
For the sake of comparison in the bottom-left panel of this figure we also show as a dashed- 
dotted line the profile obtained by Salaris et al. (1997) for a white dwarf of 0.61 Mq, a 
mass value very close to that of this panel. Note that the central oxygen abundance in the 
0.61 Mq white dwarf of Salaris et al. (1997) is somewhat higher than that of our 0.609 Mq 
white dwarf. This is mostly because we use the value of the NACRE compilation for the 
^^C(a, 7)^^0 reaction rate (Angulo et al. 1999) which is smaller than the rate of Caughlan et 
al. (1985) adopted by Salaris et al. (1997). However, another point that has to be considered 
in this comparison is the fact that the white dwarf mass in Salaris et al. (1997) is assessed 
at the first thermal pulse. Thus, for a given white dwarf mass, our progenitor stars are less 
massive than those of Salaris et al. (1997), see our Fig. [2j Because of this effect alone, our 
white dwarf model should be characterized by a higher central oxygen abundance than it 
would have resulted from a more massive progenitor — see Althaus et al. (2010). Thus, on 
the one hand we expect a lower oxygen abundance in our model because of our adopted cross 
section for the ^^C(a, 7)^^0 reaction rate, but on the other hand, we expect a higher oxygen 
abundance because of the lower initial mass of our progenitor star. The net effect is that 
the central oxygen abundance in our model results somewhat lower than that in the Salaris 
et al. (1997) model. Finally, the treatment of extra mixing episodes during core helium 
burning, which are well known to influence the final carbon oxygen stratification, leads to 
some differences in the expected composition. In our simulation, extra-mixing episodes is 
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treated as a diffusion processes (Herwig et al. 1997), while in Safaris et al. (1997) extra- 
mixing is considered as a semiconvective process. However, as shown by Straniero et al. 
(2003) both treatments give rise to a quite simiiar core chemicaf stratification, and thus 
no appreciabie difference in the central oxygen abundance should be expected from these 
treatments. 



3.7. Colors and the blue hook 

The molecular hydrogen formation at low effective temperatures also affects the evolu- 
tion of our models in the color-magnitude diagram, as shown in Fig. [HI which displays the 
run of the absolute visual magnitude My in terms of four standard colours: V — I, U — V, 
B — V, and V — R. For the evolutionary stages computed in this work the turn to the blue 
at My ~ 17 is noticeable for the V — I and V — R colors. This effect is due to the H2-H2 
collision-induced absorption over the infrared spectral regions, which forces stellar flux to 
emerge at shorter wavelenghts. Note that in this diagram, all our sequences are expected to 
become markedly blue at l^ — J and V — R. In the V — I color index, the turn-off point occurs 
between My = 16.5 and 17.2, depending on the stellar mass value. This range corresponds 
to luminosities between log(L/L0) = —4.6 and —5.0, and cooling ages between 11.9 and 15 
Gyr for our sequences with carbon oxygen phase separation. 

On the other hand, U — V and B — V colors are sensitive to the Lya broadening by 
H-H2 collisions, which reduces the emergent flux at ultraviolet and blue regions increasing 
the reddening of these colors in the models cooler than about Tgg = 5000 K { U - V > 1.0 
andB -V > 0.8). 



4. Conclusions 

We have computed a set of cooling sequences for hydrogen-rich white dwarfs, which are 
appropriate for precision white dwarf cosmo chronology. Our evolutionary sequences have 
been self-consistently evolved from the zero age main sequence, through the core hydrogen 
and helium burning evolutionary phases to the thermally pulsing asymptotic giant branch 
and, ultimately, to the white dwarf stage. This has been done for white dwarf progenitors 
with two different metallicities. For the first set of evolutionary sequences we have adopted 
solar metallicity. This allows us to obtain accurate ages for white dwarfs in the local Galactic 
disk. The second set of cooling sequences corresponds to a metallicity typical of most Galactic 
globular clusters, Z = 0.001, thus allowing to obtain accurate ages for metal-poor stellar 
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systems. To the best of our knowledge, this is the first set of self-consistent evolutionary 
sequences covering different initial masses and metallicities. 

Our main findings can be summarized as follows. First, our evolutionary sequences cor- 
rectly reproduce the observed initial-to-final mass relationship of white dwarfs. Specifically, 
our evolutionary calculations are in excellent agreement with the recent results of Salaris et 
al. (2009) for white dwarfs with solar metallicity progenitors. Second, we corroborate the 
importance of convective coupling at low luminosity in the cooling of white dwarfs, as origi- 
nally suggested by Fontaine et al. (2001). Third, we demonstrate the importance of residual 
hydrogen burning in white dwarfs resulting from low-metallicity progenitors. Fourth, we 
confirm as well the importance of carbon-oxygen phase separation upon crystallization, in 
good qualitative agreement with the results of Garcia-Berro et al. (1988a; 1988b), Segretain 
et al. (1994) and Salaris et al. (1997; 2000). Although the computed delays are smaller than 
those previously estimated by Segretain et al. (1994), they are larger than those computed by 
Salaris et al. (2000), and are by no means negligible if precision white dwarf cosmochronol- 
ogy is to be done. However, we would like to mention that these delays depend crucially 
on the previous evolutionary history of white dwarf progenitors and, particularly, on the 
rate of the ^^C(q;,7)^^0 nuclear reaction, as well as on the adopted treatment for convec- 
tive mixing. Additionally, since our evolutionary sequences rely on state-of-the-art non-gray 
model atmospheres, they reproduce the well-known blue hook of very old hydrogen-rich 
white dwarfs caused by H2-H2 collision-induced absorption (Hansen et al. 1999). Finally, 
we show the impact of Lya quasi-molecular opacity on the evolution of cool white dwarfs in 
the color-magnitude diagram. 

We would like to emphasize that our full treatment of the entire evolutionary history of 
white dwarfs has allowed us to obtain consistent white dwarf initial configurations. In par- 
ticular, the mass of the hydrogen-rich envelope and of the helium buffer were obtained from 
evolutionary calculations, instead of using typical values and artificial initial white dwarf 
models. This has implications for the cooling rates of old white dwarfs, as the thicknesses 
of these outer layers control the cooling speed of such white dwarfs. Another important 
issue which we would like to mention is that our calculations also yield self-consistent in- 
terior chemical profiles. This also has relevance for the cooling of white dwarfs, as the 
release of latent heat and gravitational energy due to carbon-oxygen phase separation upon 
crystallization crucially depend on the previous evolutionary history of white dwarfs. Also, 
the chemical stratification of white dwarf progenitors is important for correctly computing 
the specific heat of white dwarf interiors. We would like to stress as well that the evo- 
lutionary tracks of cooling white dwarfs presented here has been computed with the most 
accurate physical inputs and with a high degree of detail and realism. In particular, our 
calculations include nuclear burning at the very early phases of white dwarf evolution — 



-21 - 



which is important to determine the final thickness of the hydrogen-rich envelope — dif- 
fusion and gravitational settling — which are important to shape the profiles of the outer 
layers — accurate neutrino emission rates — which control the cooling at high luminosi- 
ties — crystallization and phase separation of carbon and oxygen — which dominate the 
cooling times at low luminosities — a very detailed equation of state — which is impor- 
tant in all the evolutionary phases — and improved non-gray model atmospheres — which 
allow for a precise determination of white dwarf colors and outer boundary conditions for 
the evolving models. Finally, we would like to remark that these evolutionary sequences are 
important as well for the calculation of self-consistent models of pulsating DA white dwarfs. 
All in all, we consider that our cooling sequences represent a landmark in the calculation 
of white dwarf evolutionary tracks, and they open an new era to precision white dwarf cos- 
mochronology. Detailed tabulations of our evolutionary sequences are available at our Web 
site |http : //www . f caglp . unlp . edu . ar/evolgroup( 
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Table 1: Initial and final stellar mass (in solar units), total masses of H and He left in the 
white dwarf (in solar units), main-sequence lifetimes (in Gyr), mass lost on the AGB and 
number of thermal pulses on the AGB for the two metallicities studied here. 
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-3.441 


-1.567 


11.885 


-0.016 
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1.00 
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-3.586 


-1.589 


9.040 


-0.104 
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0.553 


-3.577 


-1.635 


6.406 
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1.25 
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Note. — The mass of the hydrogen and hehuni contents is given at the point of maximum effective 
temperature at the beginning of the white dwarf coohng branch 



Table 2: Cooling ages when carbon-oxygen phase separation is neglected and the accumu- 
lated time delays introduced by chemical fractionation at crystallization for the evolutionary 
sequences with progenitors with Z = 0.01. 
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Table 3: Same as Table 2 but for Z = 0.001. 
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Fig. 1. — Hertzsprung-Russell diagram of our evolutionary sequences for Z = 0.01. From 
bottom to top: evolution of the 1.0 M©, 1.5 Mq, 1.75 M©, 2.0 Mq, 2.25 Mq, 2.5 Mq, 3.0 Mq, 
3.5 Mq, 4.0 Mq and 5.0 Mq model stars. Evolutionary tracks are shown from the ZAMS to 
advanced stages of white dwarf evolution. Note that the less massive sequence experiences 
a hydrogen subflash before entering its final cooling track. 
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Fig. 2. — Theoretical initial-to-final mass relationship — thick solid line — and mass of 
the degenerate core at the beginning of the first thermal pulse — thick dot-dashed-dashed 
line — obtained in this work, both for the case in which the solar composition is adopted. 
The initial-to-final mass relationships of Salaris et al. (1997) — short dashed line — of 
Dominguez et al. (1999) — long dashed hne — and of Weiss & Ferguson (2009) — dot-dot- 
dashed line — are also shown. The observational initial-to-final mass relationship of Catalan 
et al. (2008a) and Salaris et al. (2009) are the dot-dashed and solid lines, respectively. See 
the online version of the journal for a color version of this figure and the main text for 
additional details. 
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Fig. 3. — Time dependence of the different luminosity contributions for our 0.609 M© white 
dwarf sequence when carbon-oxygen phase separation is included. We show the photon 
luminosity, Lsur (solid line), the luminosities due to nuclear reactions — proton-proton chains, 
Lpp (long dashed line), CNO bicycle, Lcno (dot-dashed line), helium burning, L^e (dot- 
dot-dashed line) — the neutrino losses, Lncu (dotted line) and the rate of gravothermal 
(compression plus thermal) energy release, Lg (dashed line). Time is expressed in years 
counted from the moment at which the remnant reaches logTgfj = 4.87 at high luminosity. 
The various physical processes occuring as the white dwarf cools down are indicated in the 
figure. The progenitor corresponds to a solar-metallicity 2.0 Mq star. 
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Fig. 4. — Temporal evolution of the hydrogen content Mh (in solar masses) and the ratio 
of hydrogen (proton-proton chains and CNO bicycle) nuclear burning to surface luminosity, 
solid and dashed lines, respectively. Thick (thin) lines correspond to progenitors with Z = 
0.01 {Z = 0.001). Note that the hydrogen content left in the white dwarf, and therefore 
the nuclear energy output, are strongly dependent on the metal content of the progenitor 
stars. Although a large fraction of the hydrogen content is burnt before the remnant reaches 
the terminal cooling track at young ages, note that in the case of low metallicity, residual 
burning during the white dwarf stage reduces the hydrogen content considerably. The mass 
of the white dwarf corresponding to Z = 0.01 {Z = 0.001) is 0.609Mq (O-SOSM©). 
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Fig. 5. — Abundance by mass of H, ^He, ^^C, ^^C, ^^N, ^^O as a function of the outer 
mass fraction log(l — Mr/M^,) for our 0.609 M© white dwarf sequence at various selected 
evolutionary stages. The upper-left panel corresponds to the start of the cooling branch 
(logTeff = 4.87 at high luminosity). The upper-right panel shows the chemical profiles 
after some diffusion has already taken place in the envelope. The bottom-left panel shows 
the situation at the domain of the pulsating DA white dwarfs. Finally, the bottom-right 
panel shows the chemical abundances after the onset of cyrstallization. Luminosity and 
effective temperatures [log L/Lq, logTeff) are specified for each stage. The metallicity of 
the progenitor star is Z = 0.01. 
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Fig. 6. — Some evolutionary properties corresponding to our 0.878 M© and 0.525 Mq white 
dwarf sequences with carbon-oxygen phase separation resulting from Z = 0.01 progenitors. 
We show in terms of the surface luminosity the run of the cooling times (dashed lines), of 
the central temperature (dot-dashed lines), both read on the left-hand-side scale, and the 
evolution of the growth of the crystallized core (dotted lines) as given by the outer mass 
fraction on the right-hand-side scale. The gray area marks the occurrence of convective 
coupling. The energy released due to convective coupling and the energy resulting from the 
crystallization process markedly impact the cooling curves, see text for details. 
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Fig. 7. — Cooling curves at advanced stages in the white dwarf evolution for our sequences 
of masses 0.525 Mq (upper-left panel), 0.570 M© (upper-right panel), O.GOQM© (bottom-left 
panel), and 0.877Mq (bottom-right panel). Solid lines correspond to the case in which both 
latent heat and carbon-oxygen phase separation are considered, while dashed lines corre- 
spond to the situation when carbon-oxygen phase separation is neglected. The metallicity 
of progenitor stars is Z = 0.01. 
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Fig. 8. — Oxygen abundance (by mass) profiles for four selected white dwarfs sequences with 
stellar masses (in solar masses) as labeled in each panel. Dashed lines show the abundance 
distribution shortly after the AGB stage. Solid lines show the chemical profile after Rayleigh- 
Taylor rehomogenization has occurred in the core and dot-dot-dashed lines show the profiles 
after carbon-oxygen phase separation upon crystallization has finished. In the interests of 
comparison, in the bottom-left panel we also show with a dash-dotted line the oxygen profile 
of Safaris et al. (1997) after the rehomogenization phase of a white dwarf with mass 0.609 Mq. 
Note that in all cases, the initial chemical distribution in the core has been markedly altered 
by the crystallization process. The metallicity of the progenitor stars is Z = 0.01. 
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Fig. 9. — Absolute visual magnitude My in terms of the color index V ~ I (upper- left 
panel), U — V (upper-right panel), B — V (bottom-left panel), and V — R (bottom-right 
panel), for the complete evolutionary tracks of our sequences with Z = 0.01, from top to 
bottom: 0.525 Mq, 0.570 Mq, 0.609 Mq, 0.659 Mq, 0.767 M©, and 0.878 M©. 



